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ABSTRACT 

The  interaction  of  two  impinging  plane  turbulent  jets  has  been 
investigated  for  the  purpose  of  understanding  the  performance  of  propor- 
tional fluid  amplifiers  and  momentum-exchange  devices.   Velocity  and 
turbulence  profiles  were  measured  across  and  at  various  distances  along 
the  combined  jet  axis  through  the  use  of  a  constant  temperature  hot-wire 
anemometer.   Experimentally  determined  jet  deflection  angles  were  found 
to  agree  quite  well  with  those  predicted  in  a  previous  work  by  the  free- 
streamline  theory.   The  results,  which  are  presented  in  terms  of  normalized 
parameters,  show  that  the  jet  interaction  may  be  divided  into  three  re- 
gions:  mixing;  transition;  and  fully  developed.   A  vena-contracta  ef- 
fect was  found  to  exist  in  the  mixing  region.   A  discussion  of  each  of 
the  three  regions  is  presented  in  detail  along  with  generalizations  for 
similar  geometries. 
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I.   INTRODUCTION 

In  recent  years ,  the  need  for  a  thorough  understanding  of  the  inter- 
action of  two  or  three  jets  has  become  quite  evident  because  of  the  in- 
creased application  of  jet  interaction  in  fluid  power  technology.   There 
are  several  important  applications  which  may  be  cited  herein.   In  the 
field  of  fuel  mixing,  such  as  in  jet  or  spacecraft  engines,  the  intricate 
blending  of  two  or  more  fuel  jets  is  a  crucial  factor  in  determining  the 
engine  performance.  Another  example  is  found  in  the  area  of  fluidics. 
The  fluidic  devices  operate  on  a  principle  that  one  or  more  low  energy 
control  jets  can  impart  sufficient  deflection  to  a  power  jet  and,  thus, 
cause  a  significant  change  between  the  power  levels  of  the  receivers. 

The  proportional  fluid  amplifier  makes  use  of  the  most  direct  appli- 
cation of  the  interaction  of  three  jets.  Another  fluidic  element  or 
logic  device  called  "AND"  gate  operates  through  the  interaction  of  two 
normally  impinging  jets.   A  simplified  diagram  of  a  proportional  ampli- 
fier is  shown  in  Fig.  1.   The  deflection  of  the  main  jet  is  achieved  by 
an  imbalance  in  the  velocities  or  power  levels  of  the  two  symmetrical 
control  jets.   Whenever  a  differential  pressure  exists  between  the  control 
ports,  the  power  jet  is  deflected  by  the  momentum  imparted  by  the  control 
jets.   This  leads  to  a  large  or  amplified  pressure  signal  between  the 
receiver  or  load  ports.   The  interaction  region  of  the  amplifier  is 
vented  in  order  to  avoid  wall  attachment  due  to  the  Coanda  effect.   Thus, 
the  only  external  influence  on  the  main  jet  is  due  to  control  flows. 

Although  discovered  in  1959,  very  little  is  actually  known  about  the 
fluid  dynamics  of  proportional  amplifiers.   A  number  of  studies  have  been 
conducted  on  the  overall  performance  of  such  amplifiers,  but  the  general 
approach  in  these  investigations  has  been  restricted  to  a  "black  box" 
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treatment  rather  than  to  a  detailed  understanding  of  the  internal  mechanism 
of  the  jet  interaction*   The  realization  for  the  need  for  a  better  under- 
standing  of  the  fluid  dynamics  of  proportional  or  momentum -exchange  de- 
vices has  resulted,  in  the  last  few  years,  in  the  publication  of  several 
works  (Refs.  1-4). 

Some  basic  research  on  jet  interaction  has  been  undertaken  by  Professor 
T.  Sarpkaya  in  the  Department  of  Mechanical  Engineering  of  the  Naval  Post- 
graduate School.   The  present  research  is  part  of  the  investigation 
sponsored  by  the  U.  S.  Army  Research  Office,  Durham,  North  Carolina. 

The  present  thesis  is  a  continuation  of  a  project  undertaken  by  N. 
Gungor  (Ref.  5)  to  examine  the  interaction  of  two  impinging  jets.   The 
evolution  of  the  mixing  region  and  the  jet  deflection  angles  were  studied 
through  the  use  of  a  constant  temperature  hot-wire  anemometer.   The  veloc- 
ity and  turbulence  intensity  profiles  in  the  combined  jet  were  obtained 
for  various  control  and  power  jet  velocity  combinations. 
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II.   jBCPERIMENTAL  EQUIPMENT 

The  experimental  apparatus  together  with  the  instrumentation  used  is 
shown  in  Fig.  2.   The  basic  components  of  the  equipment  were:  (a)  jet 
assembly;  (b)  supply  system;  (c)  hot-wire  anemometer  system;  and  (d)  the 
hot-wire  calibrator. 

A.    JET  ASSEMBLY 

The  plexiglass  jet  assembly  is  shown  in  Fig.  3.   In  order  to  model 
part  of  a  proportional  amplifier,  the  axes  of  the  two  jets  were  made  to 
intersect  at  90  degrees.   The  jet-assembly  design  was  essentially  similar 
to  that  used  by  Gungor  (5)  with  the  exception  of  some  additional  modifi- 
cations to  increase  the  accuracy  and  ease  of  measurements. 

The  basic  construction  consisted  of  two  10  inch  long,  1  inch  deep, 
and  1/4  inch  wide  channels  sandwiched  between  two  1/2  inch  plexiglass 
plates.   The  control  flow  channel  had  a  1/4  inch  setback  as  shown  in  Fig. 
6.   The  outer  wall  of  each  channel  was  cut  to  a  sharp  point  (1/32  inch) 
to  prevent  attachment  of  the  jet  due  to  the  Coanda  effect.  All  sections 
were  firmly  assembled  with  dowels,  and  critical  dimensions  set  with  gage 
blocks  to  an  accuracy  of  0.001  inch.   The  channel  surfaces  were  carefully 
hand-polished  with  rouge. 

The  supply  lines  for  each  channel  were  3/4  inch  diameter  Tygon  tubing 
connected  to  special  inlet  devices  to  provide  for  the  transition  from  cir- 
cular to  rectangular  flow  channels.   Wire  screens  were  attached  between 
the  tubing  and  the  inlet  device  by  a  flange  connection.   These  screens 
served  a  dual  purposes  (1)  Eliminate  any  swirl  that  may  have  been  im- 
parted to  the  flow  in  the  supply  tubing;  and  (2)  Filter  any  foreign  matter 
out  of  the  flow. 
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The  top  and  bottom  plates  had  two  slots  cut  out  to  enable  the  hot- 
wire anemometer  mechanism  to  traverse  along  lines  parallel  to  the  combined 
jet  axis.   A  telescope  and  slide  mechanism  were  attached  to  the  bottom 
plate.   Degree  marks  were  scribed  (away  from  the  interaction  region)  into 
the  bottom  plate,  and  a  pointer  was  attached  to  the  telescope  slide. 

B.  SUPPLY  SYSTEM 

Air  was  supplied  to  the  system  from  a  50  Hp.  compressor  at  approxi- 
mately 200  psig.   The  air  was  first  passed  through  a  microfilter  to  remove 
oil,  dust,  and  other  impurities  from  the  flow  and  then  through  the  pri- 
mary pressure  regulator  (maximum  input  400  psig  and  maximum  output  of  125 
psig.).   Both  the  main  and  control  jets  were  supplied  air  independently 
from  this  primary  pressure  regulator. 

Each  jet  supply  line  had  a  series  of  pressure  regulators  (10  to  130 
psig.)  providing  an  accurately  regulated  flow  into  a  Fisher-Porter  rota- 
meter.  Both  the  main  and  control  rotameters  were  identical  having  a  maxi- 
mum capacity  of  14.6  scfm  at  14.7  psia  and  70°F  with  an  FP-3/4-27-G-10 
tube  and  a  3/4-GNSFGT-54  float.   Since  actual  flow  rates  were  not  measured 
at  14.7  psia,  each  rotameter  was  equipped  with  a  pressure  gage,  and  the 
flow  rates  were  obtained  by: 

Q  -  (rotameter  reading  (%   of  full  flow))  x  (14.6)  x  CF(P) 
where  CF(P)  -  pressure  correction  factor  obtained  from  a  chart 
provided  by  the  Fisher-Porter  company  for  the  rotameter. 

C.  HOT-WIRE  ANEMOMETER  SYSTEM 

The  Model  1050  constant  temperature  hot-wire  anemometer  manufactured 
by  the  Thermo-Systems,  Inc.  was  used  in  the  measurement  of  velocities  and 
turbulence.   The  hot-wire  probe  used  in  the  experiments  was  also  manufactured 
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by  the  Thermo-Systems,  Inc.  (wire  #  1272-T1.5).   The  wire  had  a  cold 
resistance  of  7.07  ohms  and  an  operating  resistance  of  10.61  ohms. 

A  special  probe  holder  was  constructed  and  is  shown  in  Fig.  4.   This 
holder  fit  into  slots  in  the  jet  assembly  which  enabled  the  probe  to  move 
along  lines  parallel  to  the  combined  jet  axis.  A  micrometer  mechanism  on 
the  holder  provided  transverse  movement  of  the  hot-wire  across  the  jet 
with  an  accuracy  of  0.001  inch. 

D.    HOT-WIRE  CALIBRATOR 

The  hot-wire  probe  was  calibrated  with  the  Thermo-Sy stems,  Inc.  Model 
1125  Calibrator  shown  in  Fig.  5.   The  calibrations  were  made  with  the 
high  range  calibration  orifice  having  a  diameter  of  0.150  inch.   The 
pressure  in  the  calibration  chamber  was  measured  with  water  manometers 
with  various  ranges  of  pressure  to  increase  the  accuracy  of  the  calibra- 
tion.  The  use  of  the  calibrator  will  be  discussed  in  the  next  section. 

A  schematic  drawing  of  the  entire  experimental  system  is  shown  in 
Fig.  7. 
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III.   EXPERIMENTAL  PROCEDURE 

The  experimental  procedure  consisted  of  the  following  six  main  steps: 
(a)  Calibration  of  the  hot-wire  probe;  (b)  Determination  of  the  approximate 
jet  deflection  angles;  (c)  Adjusting  of  the  hot-wire  anemometer  instru- 
mentation; (d)  Setting  of  the  selected  flow  rates  for  each  jet;  (e) 
Measurement  of  local  mean  velocities  and  turbulence  intensities;  and  (f) 
Data  reduction. 

A.    CALIBRATION 

The  calibration  consisted  of  the  establishment  of  the  exact  relation- 
ship  between  the  anemometer  voltage  outputs  and  the  corresponding  veloc- 
ities.  The  Model  1125  Calibrator  was  used  for  this  purpose. 

Firsts,  the  hot-wire  probe  was  placed  about  1/16  inch  above  the  ori- 
fice of  the  calibrator  as  recommended  by  the  manufacturer  (Ref.  6).  Air 
was  then  let  through  the  orifice s   and  the  pressure  in  the  orifice  chamber 
was  measured  by  a  manometer.   From  these  pressure  readings 9  the  velocity 
through  the  orifice  was  calculated  by; 


where  P  is  the  chamber  pressure  in  pgig.9  and  />  is  the  air  density. 
For  each  velocity  the  corresponding  voltage  output  of  the  hot-wire  anemo- 
meter was  recorded.   Finallys  a  calibration  curve  of  velocity  versus  volts 
was  plotted  as  shown  in  Fig.  8. 

The  calibration  curve  in  this  form  is  not  very  useful  since  it  cannot 
be  used  easily  for  a  large  number  of  measurements.   In  addition  the  de- 
rivatives of  the  calibration  curve,  which  are  needed  for  the  turbulence 
intensity  calculations s  cannot  be  obtained  from  a  calibration  curve  of 


this  type.  A  polynomial  approximation  is  obviously  necessary ,  especially 
for  computerized  data  reduction„ 

The  hot-wire  anemometer  characteristics  follow  King's  Law  (Ref.  7) 
which  is  essentially  a  power  law  relation  between  voltage  and  velocity. 
A  plot  of  log  (velocity)  versus  log  (voltage)  data  is  shown  in  Fig.  9  and 
is  seen  to  be  very  near  to  a  straight  line.  A  second  degree  polynomial 
approximation  of  this  curve  was  found  to  fit  almost  exactly  to  the  logari- 
thmic form  of  the  data  points.   The  derivation  of  this  equation  is  given 
in  Appendix  A.   The  polynomial  approximation  to  the  calibration-curve  is 
shown  in  Fig.  10  together  with  the  actual  data  points  for  comparison. 

The  calibration  curve  described  above  was  checked  frequently  during 
the  course  of  the  experimentation.   It  was  found  that  the  wire  performance 
remained  essentially  the  same. 

B.   APPROXIMATE  JET  DEFLECTION  ANGLE 

In  order  to  properly  place  the  hot-wire  probe  near  the  center  of  the 
jet,  the  approximate  jet  deflection  angle,  oC  ,  had  to  be  determined. 
Two  approximate  methods  were  used  by  Gungor  (5),  but  they  were  found  to  be 
rather  time  consuming  and  tedious.   Since  the  accuracy  required  was  only 
+  2  degrees,  a  simple  needle  and  thread  method  was  used  to  determine  oc     . 

The  apparatus  consisted  of  a  plastic  holder  with  a  long  needle  at- 
tached. At  the  end  of  the  needle  an  eight  inch  piece  of  yarn  was  care- 
fully glued  so  as  to  have  no  preferred  directional  effects  due  to  the 
method  of  attachment.   The  needle  end  was  then  placed  at  the  center  of 
the  nozzle  (Point  M  in  Fig.  11).   For  each  main  and  control  flow  combina- 
tion, the  approximate  jet  angle  could  easily  be  read  from  the  angle  the 
thread  made  with  the  axis  of  the  main  jet. 
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C.  ADJUSTING  OF  ANEMOMETER  INSTRUMENTS 

Prior  to  each  experiment,  the  anemometer  was  turned  on  and  allowed 
to  warm  up  for  about  one-half  hour.   The  cold  resistance  of  the  wire  was 
recorded  and  checked  to  make  sure  that  it  had  not  deviated  from  previous 
values.   Next  the  operating  resistance  was  set  into  the  device  in  accord" 
ance  with  the  operating  instructions  (Ref.  6).   The  system  was  then  switched 
to  the  "RUN"0  operation  mode„ 

The  fluctuating  signal  was  monitored  by  an  oscilloscope  throughout 
the  experiment.   The  internal  feedback  stability  was  adjusted  by  the  "TRIM" 
dial  on  the  anemometer  to  render  the  system  stable  within  the  range  of 
velocities  used. 

D.  SETTING  OF  FLOW  RATES 

The  channel  flow  rates  were  measured  by  the  main  and  control  rota- 
meters and  adjusted  by  two  pressure  regulators.   The  main  jet  was  main- 
tained at  an  average  velocity  of  87.5  ft/sec  which  corresponded  to  a 
rotameter  setting  of  60%  of  full  flow  and  a  Reynolds  number  of  11,140. 
The  control  jet  flow  rate  was  varied  in  the  experiment  to  obtain  differ- 
ent deflection  angles.   The  control  flow  rates  used  and  the  corresponding 
velocities  and  Reynolds  numbers  are  tabulated  below. 
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Control  Flow  Rate 
(7.  full  flow) 

scfm 

U 
c 

ft/sec 

Reynolds 

Number 

U  /U 

c   m 

60 

9.11 

87.5 

11,140 

1.00 

50 

7.59 

72.9 

9,287 

0.834 

40 

6.01 

57.7 

7,340 

0.66 

30 

4.47 

42.9 

5,465 

0.491 

20 

2.95 

28.3 

3,600 

0.323 

10 

1.46 

14.02 

1,786 

0.1605 

NOTE;   Reynolds  numbers  are  based  on  the 
width  of  the  main  channel. 

E.   MEASUREMENT  OF  VELOCITY  AND  TURBULENCE  INTENSITIES 

Velocity  and  turbulence-intensity  measurements  were  made  along  lines 
perpendicular  to  the  jet  axis  at  nine  different  distances  from  the  inter- 
section  of  the  inner  walls  of  the  jet  channels.  These  nine  stations  were 
chosen  to  give  both  an  accurate  picture  of  the  mixing  region  and  the  fully 
developed  region.  The  following  station  distances  were  used;  0.15  inch; 
0.30  inch;  0.45  inch;  0.60  inch;  0.80  inch;  1.00  inch;  1.50  inchs;  2.50 
inches;  and  4.00  inches. 

The  probe  was  positioned  using  the  telescope  assembly  discussed  pre- 
viously.  The  telescope  was  first  set  to  the  approximate  jet  deflection 
angle.   The  probe  was  then  centered  at  a  particular  station  with  the  use 
of  the  telescope  crosshairs. 

Once  the  probe  was  positioned,  the  traversing  limits  were  determined 
to  give  an  accurate  picture  of  the  jet  profile.   The  hot-wire  probe  was 
then  moved  at  suitable  increments  along  the  traversing  line,  and  both  the 
wire  voltages  and  the  R.M.S.  voltages  were  recorded. 

A  schematic  drawing  of  the  traversing  line  together  with  the  appropri- 
ate nomenclature  is  shown  in  Fig.  11. 


F.    DATA  REDUCTION 

All  data  was  reduced  and  plotted  using  the  XBM°360  computer  <>   The 
main  data  reduction  program  is  presented  in  Computer  Program  B.  Veloc' 
ities  and  turbulence  intensities  were  calculated  using  the  polynomial 
calibration-curve  approximation  previously  described. 


The  turbulence  intensity  was  calculated  by: 

dE 


U   U 


where 


fe   is  the  R.M,S.  voltage  reading  and  dU  is  the  derivative  of 

dE 
the  calibration  curve  evaluated  at  the  value  of  the  local  mean  velocity 


U„   The  derivation  of  the  above  equation  is  given  in 


ix 


IV.   DISCUSSION  OF  RESULTS 

The  main  purpose  of  this  investigation  was  to  gain  a  better  under- 
standing of  the  mixing  of  two  turbulent  jets.   This  interaction  was  found 
to  be  dominated  by  three  regions:   (1)  the  mixing  region;  (2)  the  transi- 
tion region;  and  (3)  the  fully  developed  region.   The  important  character- 
istics of  each  of  these  regions  will  be  discussed  in  terms  of  the  velocity 
and  turbulence  intensity  measurements.  General  similarity  in  jet  inter- 
action for  similar  geometries  will  then  be  discussed  through  a  comparison 
of  the  data  presented  herein  with  those  obtained  by  Gungor  (5). 

The  most  important  characteristic  in  any  jet  interaction  study  is 
obviously  the  jet  deflection  angle.   Figure  12  shows  the  comparison  of 
experimental  jet  deflection  angles  with  those  predicted  by  the  free- 
streamline  theory  (5).  Also  plotted  on  this  graph  are  the  jet  deflection 
angles  based  on  a  simple  momentum  equation: 

1   /  U  \  2 
©C=  tan"     jp\ 

\    m  I 

It  is  apparent  that  the  predictions  of  the  free-streamline  analysis 
and  the  experimental  values  agree  quite  satisfactorily.  Deflection 
angles  based  on  the  momentum  equation  do  not  agree  well  at  all.   This  is 
to  be  expected  since  the  equation  does  not  take  into  account  the  control 
channel  setback. 

Gungor  (5)  also  measured  the  deflection  angles  through  the  use  of  a 
model  half  the  size  of  that  used  in  the  present  investigation,  i.e., 
W  =  1/8  inch,  aspect  ration  *  8.   The  table  below  gives  a  comparison  of 
the  two  sets  of  deflection  angles.   Both  sets  of  data  have  similar  geome- 
tries and  comparable  Reynolds  numbers  with  the  exception  of  the  aspect 
ratios. 
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Flow  Ratio 
c  m 

6/6 

5/6 

4/6 

3/6 

2/6 

1/6 

W  ■  1/4  in. 

W  =  1/8  in. 
(Gungor) 

35 
35 

30 
30 

25 
25 

18.5 
18 

11 
10 

2 
2.5 

It  is  obvious  that  the  deflection  angles  are  almost  identical  al- 
though the  actual  magnitudes  of  the  velocities  were  quite  different  in 
each  case.  Apparently,  the  deflection  angle  depends  only  on  the  velocity 
ratio  rather  than  the  actual  magnitudes  of  the  velocities. 

The  three  regions  of  the  jet  interaction  can  be  seen  clearly  in  Figs. 
13-21.   In  these  graphs  the  velocity  and  turbulence  intensity  profiles 
for  equal  main  and  control  jet  velocities  are  plotted  for  various  dis- 
tances along  the  jet  axis.   The  mixing  region  is  near  the  nozzle  exit, 
and  the  two  flows  remain  essentially  separate  with  two  velocity  peaks. 
At  larger  distances  from  the  exit,  the  two  jets  mix  together  until  a 
single  symmetrical  jet  emerges  at  a  distance  of  about  6W  as  shown  in  Fig. 
19.   This  is  the  fully -developed  region.   The  area  between  the  mixing  and 
the  fully-developed  regions  is  the  transition  region. 

Also  shown  in  this  set  of  graphs  are  the  turbulence  intensity  pro- 
files.  This  combination  of  velocity  and  turbulence  profiles  shows  graphi- 
cally the  two  turbulent  mixing  regions  at  the  outer  edges  of  the  jet. 
These  mixing  regions  are  caused  by  the  free  shear  layers.   A  third  turbu- 
lent mixing  region  exists  in  the  central  core  between  the  two  jets  and 


will  be  discussed  later. 
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A  complete  set  of  data,  for  nine  stations  along  the  jet  axis,  is 

presented  for  U  /U   =1.00.   Only  two  representative  velocity  and  turbu- 
c  m 

lence  intensity  profile  combinations  are  presented  for  each  of  the  remain- 
ing velocity  ratios  (Figs.  22-31).   This  is  because  all  of  the  control 
flows  exhibited  similar  characteristics.   Turbulence  intensities  for 

U  /U  =  0.1605  have  not  been  presented  because  of  the  inaccuracies  in  the 
c  m 

derivatives  of  the  calibration  curve  for  low  velocities. 

The  mixing-region  normalized-velocity  profiles,  for  each  flow  ratio, 
are  presented  in  Figs.  32-37.   From  these  graphs  the  actual  mixing  pro- 
cess of  the  jets  can  be  seen  quite  clearly.   Close  to  the  channel  exit, 
the  two  velocity  maximums  are  far  apart  with  a  large  velocity  minimum  or 
"trough"  between  them.   This  preservation,  for  a  short  distance,  of  the 
identity  of  the  main  and  control  jets  has  been  attributed  by  Douglas  and 
Neeve  (2)  to  a  "brick  wall  effect".   Further  down  the  jet  axis,  the  peaks 
move  closer  together  as  the  momentum  transfer  across  the  turbulent  shear 
layer  decreases.  Attemps  were  made  to  find  some  kind  of  similarity  vari- 
able for  the  mixing-region  velocity  profiles  but  to  no  avail.   Similarity 
will  probably  be  obtained  by  analyzing  each  velocity  mixing  region  sep- 
arately as  mentioned  in  Pai  (Ref.  8). 

In  the  mixing  region  the  turbulence  intensities  are  higher  on  the 
control  jet  side  than  on  the  main  jet  side.   This  higher  turbulence  in- 
tensity is  caused  by  the  low  local  mean  velocities  and  not  by  the  R.M.S. 
value  of  turbulence.   With  decreasing  control  flow,  the  difference  be- 
tween the  control  and  main  jet  turbulence  intensity  peaks  rise  which  sup- 
ports the  above  reasoning.   This  same  trend  is  also  verified  by  Gungor 
(5). 
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As  previously  mentioned,  there  are  three  turbulent  mixing  regions. 
The  graphs  cited  above  presented  only  the  turbulence  intensities  in  the 
exterior  regions  of  the  velocity  profiles.   Figures  38-40  show,  on  an 
extended  scale,  the  internal  turbulence  profiles  for  three  velocity 
ratios.   The  intensity  of  internal  turbulence  decreases  with  increasing 
distance  from  the  channels,  as  would  be  expected.   By  comparing  Figs. 
38  and  40,  it  is  seen  that  the  greater  the  velocity  difference  between 
the  two  jets,  the  greater  the  internal  turbulence  intensity.   This  in- 
ternal turbulence  is  a  result  of  the  mixing  of  the  two  jets  due  to  vor- 
tex motion.   The  turbulence  decreases  because  of  the  vortex  decay  in 
the  shear  layer. 

It  is  apparent  from  the  definition  of  turbulence  intensity  (see 
Appendix  B)  that  it  also  represents  the  relative  noise  in  the  jet. 
Moynihan  (Ref.  3)  took  the  total  pressure  ratio  as  a  qualitative  measure 
of  noise  in  a  jet  interaction.   In  graphs  similar  to  those  presented 
here,  he  showed  three  distinct  noise  levels  in  the  jet  of  about  equal 
magnitude.   The  results  presented  herein  show  that  they  are,  in  fact, 


not  equal.  The  amount  of  turbulence  (  y  u  )  is  much  less  in  the  internal 
shear  layer  than  in  the  external  mixing  regions.  In  proportional  amplifi- 
er design,  this  fact  is  very  important  to  keep  in  mind. 

The  transition  region  has  been  defined  as  an  area  where  the  veloc- 
ity profiles  do  not  exhibit  dual  velocity  maximums  nor  are  they  symmetri- 
cal as  in  a  fully-developed  jet  profile.   A  good  example  of  this  type  of 
velocity  profile  is  shown  by  Curve  A  in  Fig.  43.   Using  crude  terminology 
to  explain  the  transition  process,  the  jet  ""pops"  from  the  mixing  region, 
through  the  transition  region,  and  into  the  fully  developed  region.   The 
length  of  this  transition  region,  for  all  control  flows,  was  less  than 
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one  channel  width.   This  region  is  short  because  of  the  characteristics 
of  the  vortex  decay  in  the  inner  shear  layer.  As  time  increases,  the 
turbulent  vortices  decay  very  quickly,  causing  a  very  short  transition 
before  a  fully  developed  jet  appears. 

In  Figs.  41-46,  similarity  profiles  are  presented  for  the  fully- 
developed  velocity  region.   These  graphs  confirm  the  findings  of  Gungor 
(5)  that  fully  developed  profiles  are  similar  when  plotted  against  x/D. 
For  convenience  the  transition  region  profiles  have  also  been  included 
in  these  plots.   Since  transition  profiles  are  not  symmetric,  they  can- 
not be  expected  to  follow  such  a  similarity  trend. 

The  effective  pivot  point  of  the  jet,  for  each  flow  ratio,  was 
examined  by  plotting  jet  maximum  points  along  the  jet  axis  similar  to 
that  shown  in  Fig.  11.   Within  the  accuracy  of  this  experiment,  the  pivot 
point  was  found  to  be  approximately  at  the  intersection  of  the  centerlines 
of  the  main  and  control  channels  (Point  M  in  Fig.  11)  for  all  velocity 
ratios.   For  this  reason,  all  of  the  deflection  angles  were  measured 
relative  to  this  point.   This  stationary  pivot  verifies  the  findings  of 
Moynihan  and  Reilly  (4). 

One  of  the  most  interesting  plots  of  the  mixing-region  velocity  data 
is  shown  in  Fig.  47.   This  graph  shows  that  the  ratio  of  the  two  velocity 
maximums  for  a  given  flow  ratio  remains  nearly  a  constant  in  the  mixing 
region.   This  fact  could  be  very  important  in  attempts  to  modify  the  free- 
streamline  theory  to  better  describe  the  mixing  region.   It  is  also  seen 
from  these  curves,  that  the  length  of  the  mixing  region  increases  with 

increasing  control  flow.   The  maximum  length  is  6W  for  U  /U  =  1.00, 

c  m 

The  dashed  lines  to  the  right  of  each  curve  represent  the  transition 
region  and  show  how  small  this  region  actually  is  (about  0.5W). 
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Figure  48  shows  the  acceleration  of  flow  and  the  presence  of  a  vena- 
contracta  in  the  mixing  region.   Here  the  ratio  of  the  largest  maximum 
velocity  in  a  particular  profile  to  the  maximum  velocity  in  the  main  jet 
channel  has  been  plotted  against  D/W,  A  maximum  occurs  at  a  distance  of 
2.5W  for  almost  all  flow  ratios.   This  may  be  interpreted  as  a  "'pinching" 
or  vena-contracta  of  the  combined  jet. 

This  vena-contracta  effect  is  caused  by  two  phenomena.   First,  the 
control  jet  "pinches"  the  main  jet  because  of  the  direction  of  the  initial 
control  jet  momentum.   Secondly,  as  the  jet  moves  further  from  the  nozzle, 
surrounding  air  is  entrained  by  the  turbulent  mixing.   Eventually,  this 
entrainment  and  mixing  causes  the  jet  to  expand  and  decrease  in  velocity. 
The  net  result  is  a  vena=contracta  type  of  velocity  profile.   The  amount 
of  pinching  of  the  main  jet  is  certainly  related  to  the  relative  momentum 
of  the  control  flow.   The  vena-contracta  is  almost  non-existant  for  low 
control  flows. 

One  of  the  primary  purposes  of  this  investigation  was  to  determine 
if  any  general  similarity  in  jet  interaction  existed  for  similar  geometries, 
In  order  to  accomplish  this,  the  results  for  W  a  1/4  inch  were  compared, 
wherever  possible,  with  the  results  of  Gungor  (5)  for  W  ■  1/8  inch. 

A  comparison  of  centerline  velocities9  for  various  distances  along 

the  combined  jet,  is  presented  in  Figs.  49=54.   The  U    /U  ratio  is 

seen  to  have  the  same  general  functional  relationship  with  D/W  even  though 

velocities  and  channel  widths  are  quite  different.   Figure  55  is  a  graph 

of  D/W  (corresponding  to  the  maximum  centerline  velocity)  versus  the  flow 

ratio.   It  appears  that  the  centerline  velocity  is  a  function  of  D/W  and 

U  /U   for  similar  geometries  and  comparable  Reynolds  numbers. 
c  m 
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in: 


V.   CONCLUSIONS 
The  following  conclusions  can  be  drawn  from  the  data  presented  here- 

(1)  Experimental  jet  deflection  angles  agree  with  the  predictions 
of  the  free-streamline  theory. 

(2)  The  jet  deflection  angle  is  independent  of  the  magnitudes 
of  the  main  and  control  jet  velocities  for  similar  geo- 
metries and  depends  only  on  velocity  ratio. 

(3)  For  two  impinging  jets,  there  are  three  turbulent  mixing 
regions  with  the  turbulence  intensity  directly  dependent 
on  velocity  differences. 

(4)  Fully-developed  jet  profiles  for  each  flow  ratio  are 
similar  when  plotted  against  x/D. 

(5)  The  effective  pivot  point  of  the  deflected  jet  is  the  same 
for  all  flow  ratios,  i.e.,  at  the  intersection  of  the 
channel  axes. 

(6)  In  the  mixing  region,  the  ratio  of  the  two  velocity  maxi- 
mums is  nearly  constant  for  each  flow  ratio. 

(7)  The  length  of  the  mixing  region  increases  with  increasing 
control  flow. 

(8)  A  vena-contracta  effect  is  present  in  the  combined  jet  at 
distances  of  about  2.5W  from  the  channel  exits. 

(9)  The  centerline  velocity  in  the  combined  jet  is  a  function 
of  only  D/W  and  velocity  ratio  for  similar  geometries  and 
comparable  Reynolds  numbers. 
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VI.   RECOMMENDATIONS  FOR  FURTHER  WORK 

The  investigation  described  herein  should  be  extended  to  the  study 
of  the  following  areas ; 

(1)  The  effect  of  different  setbacks  and  aspect  ratios  on 

the  jet  deflection,  velocity ,  and  turbulence  distributions , 

(2)  Effect  of  low  Reynolds  number  or  laminar  flow  on  the  jet 
interaction. 

(3)  A  more  thorough  investigation  of  the  turbulence  in  the 
jets  using  two-wire  correlation  to  check  for  dominant 
noise  frequencies  and  self -preserving  turbulence. 

(4)  Verification  of  theoretical  predictions  for  similarity  in 
the  three  velocity  mixing  regions. 

(5)  An  extension  of  the  free -streamline  theory  approximation 
for  the  mixing  region  of  two  impinging  jets  to  include 
two  unequal  free  stream  velocities  of  constant  ratio. 
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APPENDIX  A 

LOGARITHMIC  CALIBRATION  CURVE  APPROXIMATION 
The  basic  relationship  between  the  voltage  and  velocity  for  a  hot- 
wire anemometer  is  given  by  King's  Law; 


(1)    AE  =  K  Un 


where 


AE  -  differential  voltage 
U~  =  local  mean  velocity 
K,n  constants 


taking  the  logarithm  of 

(2)   ln(4E)  =  ln(K)  +  nln(U) 
letting 

y  =  ln(AE) 

V  -  ln(IT) 


and  approximating  (2)  with  a  polynomial  one  has 
(3)  V  =  ax  +  a2y  +  a3y  +„...+  am+]ym 

Various  values  of  m  were  tried  using  a  Chebyschev  Curve-Fit  as  shown 

in  Computer  Program  A.   A  second  degree  (m  -  2)  polynomial  approximation 

was  found  to  give  the  best  fit  to  the  actual  data  (see  Fig.  10). 

The  coefficients  a   for  the  polynomial  were  found  to  be: 
m 

a 

m 


a.  -  3.494994 
1 


a.  ^  3.008039 


a  »  -0. 003591135 
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DERIVATION  OF  TURBULENCE-INTENSITY  EQUATION 
For  a  turbulent  flow,  the  instantaneous  value  of  the  anemometer  out- 
put is  given  by; 

(1)   E  -  ¥+Ve 


where 


E  ■  voltage  output  corresponding  to  the  mean  velocity 

at  a  given  point 


7  e  -  R.M.S.  value  of  the  voltage  output  corresponding 
to  the  R.M.S.  value  of  the  fluctuating  component 
of  velocity  (  ~\Jw^)  . 

Thus  from  equation  (1)  we  can  write  the  equation  for  the  instantaneous 

velocity 

(2)  U  -  TJ  +"J"~2 

and 

(3)  E  -  f(U)  =  f(u  +  "yu) 

also 

(4)  dE  -  f'(U)dU 

since  the  fluctuating  component  of  velocity,   yu   ,  is  relatively  small, 
we  can  expand  equation  (3)  in  a  Taylor  Series  neglecting  higher  order 
terms, 


(5)    f(  u  -fllu  2)  -  f(TT)  +"Vu  'feu;  +  ... 

therefore 

(6)  E  -  1"  +  Ve   -  f(U)  +Yu  f ' 

_ 
by  definition  E  -  f(U)  so  that 

(7)  ifF5"  -ifu^f'CU) 
and 
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(8)  i/T1  =  VT1 

f'CO) 


the  turbulence  intensity  is  defined  as 
(9) 


T^      T^ 


0       U  f'(U) 

since  f'(U)  =  d| 

dU 

the  final  expression  for  turbulence  becomes 

<10>  izl  -  Vzl  •  s. 

U        U        dE 
Equation  (10)  is  easily  adaptable  for  computerized  data  reduction 
as  shown  in  Computer  Program  B.  The  values  of  U,  and  dll/dE  are  obtained 
from  the  polynomial  calibration  curve. 
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FIGURE  1.   PROPORTIONAL  AMPLIFIER 
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FIGURE  6.   SCHEMATIC  OF  JET  ASSEMBLY 
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M  «  Jet  Pivot  Point 
x  -  Traversing  Distance 
D  •  Station  Distance 
•c-  Deflection  Angle 


Control  Jet 
Channe 1 


FIGURE  11.   DIAGRAM  OF  TRAVERSING  LINE 
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